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A many body theory for a two-component system of spin polarized interacting fermions in a one- 
dimensional harmonic trap is developed. The model considers two different states of the same 
fermionic species and treats the dominant interactions between the two using the bosonization 
method for forward scattering. Asymptotically exact results for the one-particle matrix elements at 
zero temperature are given. Using them, occupation probabilities of oscillator states are discussed. 
Particle and momentum densities are calculated and displayed. It is demonstrated how interactions 
modify all these quantities. An asymptotic connection with Luttinger liquids is suggested. The 
relation of the coupling constant of the theory to the dipole-dipole interaction is also discussed. 
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I. INTRODUCTION 



The achievement of Bose-Einstein condensation in dilute ultracold gases (l) renewed the interest in fermionic many 
body systems and their supcrfluid properties 0-0]. Recent experimental successes in obtaining degeneracy in 
three dimensional Fermi vapors |9|JlO|] intensified the interest in confined Fermi gases. 

Using microtrap technology ]Tl|-|l4|, it will become possible in the near future to produce a neutral ultracold quantum 
■ gas of quasi one-dimensional degenerate fermions. 

In many cases, identical spin polarized fermions experience only a weak residual interaction because s-wave scattering 
is forbidden. This restriction does not hold for a two-component system of spin polarized fermions and significant 
interactions between the components are possible. For instance, the dipole-dipole interaction jll| can become relevant, 
especially in the case of polar molecules [[L6| . 

The confinement of a trapped ultracold gas can be realized by a harmonic potential, which is more realistic than 
trapping between hard walls ("open boundary conditions"). The latter system of interacting one-dimensional fermions 
constitutes a bounded Luttinger liquid, which allows an exact treatment for certain types of interactions Jl7] ^q] . 
In this article, we consider a quasi one-dimensional spin polarized Fermi gas composed of an equal number of atoms 
in two different internal states and confined by a harmonic potential. One possible realization consists of trapped 
(electron) spin polarized fermions in two different hyperfine states as discussed in || for the case of 6 Li. 
We consider the inter-component interaction between the two components. As in [p4[ , we apply the bosonization 
method known from the Luttinger model (for reviews we refer to [^5|-^|) to treat the interactions. The generalization 
to two components is analogous to the inclusion of spin fi/2 into the Luttinger model pq| . The bosonization method 
• ■ relies on fermion-boson transmutation in one spatial dimension: Physical quantities can be calculated in a bosonic 
formulation instead of the fermionic theory, and the two calculations give the same answer |29|-|3^| . 
We show how interactions modify the one-particle properties of the two-component Fermi gas. Results for the non- 



interacting Fermi gas in a one-dimensional harmonic trap were given in [ pSP'j 
The paper is organized as follows. Section II develops the theory for the two-component case. Section III applies the 
theory to the calculation of the one-particle matrix elements. In Section IV, occupation probabilities, off-diagonal 
matrix elements, and densities of particles and momenta are evaluated numerically for two different interaction models. 
Section V discusses the problem of the Fermi edge in the present case and a possible relation to the standard Luttinger 
model result is pointed out. In Section VI, the relation of the coupling constant of the theory to the dipole-dipole 
interaction is discussed. An Appendix is concerned with the bosonization procedure for bilinear forms of auxiliary 
fields in the case of two components. 



II. TWO-COMPONENT THEORY 



The two component Fermi gas of uniform mass tua is confined by the one-dimensional harmonic potential 

V(z) = \m A uJ\z\ (1) 
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with longitudinal trap frequency oj#. The unperturbed Hamiltonian in second quantization is 

oo 

H = ^ ^nC^' ( 2 ) 
n=0,cr=±l 

The index a — ±1 refers to the two components and c^ a creates a fermion of species a in the oscillator state \n). 
The one-particle energies 

hu n = fuJi(n+l/2), n = 0,l,..., (3) 

are seen to depend linearily on the quantum number n of oscillator states. This is one of the requirements for 
bosonization. In addition, exact solvability rests on the presence of an anomalous vacuum (cf. p5|-27j]), which is 
constructed by extending the linear dispersion of oscillator states to arbitrarily negative energies and then filling 
all states of negative energy. However, the anomalous vacuum has little effect for processes near the Fermi energy 
ep = hidi(N — 1/2) provided N is sufficiently large |3^| , making the treatment asymptotically exact. 
The success of the Luttinger model is based on the possibility to express forward scattering processes entirely in terms 
of the density fluctuation operators. For the two component system, these operators are 

pAp) = J2c+ +pa c q<7 . (4) 
i 

Due to the presence of the anomalous vacuum, they obey bosonic commutation relations 

lPa(-p),Pa>(q)] =P<W<W (5) 

In our case, the interaction Hamiltonian is given by a two particle interaction 

V = \ E V(m<j',pa;q<j',na) (c+ (7 ,Cq (7 ')(c+ (7 c n(T ) (6) 

m,npq,a.a' 

without "component flip", i.e., the possibility for a fermion to change its state a in the collision process is excluded. 
The case a — a! corresponds to a weak intra-component interaction Vji, while a = —oi is a relevant inter-component 
interaction V±. 

Similar to (24j| , two cases of solvable forward scattering processes can be identified: 

V = V a + V b (7) 

with 

v* = \ E v -\\ (H) Ar(-p) pAp) + \ E y aA\p\) p-A~p) pAp), (8) 

P,<7 p,(T 



K = IJ2 v "\\ (W) pAp) pAp) + \ E v ^\p\) pAp)p-Ap)- 



2^ ■»liwi^»w;ruvry g 

p,CT p. (7 

The coupling functions V a ± and Vb± are the analogues of 174^ and g-2±_ of the Luttinger model. 

Forward scattering dominates when the pair interaction is sufficiently long ranged, e.g., in the case of dipole-dipole 
interactions. In [p5| , a detailed discussion is given, how the assumed forms (0) can be related to real scattering 
potentials. 

The case of two components requires canonical transformations to mass fluctuation operators 

P(P) = -^[P+(P) + P-(P)}, (9) 



and component fluctuation operators 



a(p) = ±= 0+(p)-p_(p)], (10) 
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such that the Hamiltonian for low lying excitations separates into H = H p + H a in analogy to the spin 1/2 case E5 
of the Luttinger model. 

The transformation to new bosonic operators 



leads to canonical commutation relations 

[^m/ij^)ii/] — 3 &rn.n- (1^) 

The new label v — ±1 refers to mass and component fluctuations. 

The bosonic version of the unperturbed Hamiltonian in the N-fermion sector is ]3J4p6j : 



Ha = E 171 {dmvdmv + dmvd+v} . (13) 



m>0,u 

The complete bosonic interaction operator becomes 

i m[y o ||(m) + vV a ±(m)] jd+^d^ + d tm ,<2;t u ,} ( 14 ) 



^ 2 

m>u,y 



2 

m>0,z^ 

The total bosonic Hamiltonian is diagonalized in a standard way using the Bogoliubov transformation 

d m u = S~^fm V S — /mi; cosh £ m „ fmu sinh C^ mv , (15) 

with 



I m>0,w=±l J 



The transformation parameters Q mv are determined by the diagonalization conditions 



cosh(2C mi ,) 

In connection with the calculation of one-particle matrix elements, scaled coupling constants 



(16) 



tanh(2Cm^) = z — ' T/ , r— T/ ; r - 17) 
huji + V a \\ (m) + fl4_L (m) 

Finally, we arrive at the free bosonic Hamiltonian 

H = E m e mufmufmu + Const., (18) 

m>0,i> 

describing density wave excitations in the two-component Fermi gas. The excitation spectra are 

_ huji + V a \\{m) + vV a ±(m) 



(19) 



1 2 

ow = X sinh(2C m „), 7 m „ = sinh Cmi- (20) 

will appear. 

Usually, the intra-component scattering is negligible (Vy — > 0, the case Vj| ^ was considered in |24j] for the one- 
component system) and dominant forward scattering for the inter-component part results |3q] in 

V j.(m) = F 6X (m) = V(m) (21) 
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Then the simpler relations 

e m „ = fiwt y/l + 2vV{m), a mu = — = (22) 

2-^/1 + 2z/ K (m) 

hold for \V(m)\ < Tiut/2. 

Following |24| , we will also consider two specific interaction models: A simplified model IM1, when only one mode 
Vim) — V(l)(6 m ,i + contributes. This model preserves many of the features of the interaction in the full 

model (interaction model 2, see IM2 below), when infinitely many modes are superimposed. 
In the case of IM1, the relevant coupling constants are 

Ci„ = \ artanh ( *££L ) , <*iv = ~ sinh(2Ci,), 7 i„ - \ U± + ±a\ v - 1 J . (23) 



2 " V 1 + ^(1) / ' 2 v Sll// ' " 2 

In the case of IM2, the coupling constants decay exponentially according to 

a mu = exp(-r a m/2) a „, &0u = exp(r Q /2) a\ v , Jmu = exp(-r 7 m) 70, 7o„ = exp(r 7 ) j lv . (24) 

An essential step in the actual calculation of physical quantities is the connection between fermionic operators and 
bosonic fields. The bosonization of fermion generation and destruction operators is completely solved for the Luttinger 
model Hfflll. 

In the present case, the situation is less comfortable: Apart from the above association of mass and component 
fluctuation operators with the d operators, we can bosonize only bilinear forms of an auxiliary field following the 
prescription of Q . The auxiliary field is defined by 

00 

i; aa (v)= J2 e ilv c lr7 =i> arT {v + 2n). (25) 
i— —00 

The Appendix demonstrates that the required bosonization for a two-component Fermi gas is 

^i<j{u)4>aa{v) = G N (u - v) exp j-i (J>t(u) - 4>t( v )j } exp j-i U a (u) ~ <M«)) } , (26) 

using the two-component non-Hcrmctian bosonic field 

00 

4> a («) = -i -7= e mV (dn+ + oL- )*fc{v). (27) 



. /2n 

n—l v 



The distribution- valued prefactor Gn(u) is the same as in 40 : 



N-l 



G N (u) = J2 e- il ^ u+ir,) . (28) 



III. ONE-PARTICLE MATRIX ELEMENTS 

The above prescription allows to calculate analytically all m-particle matrix elements of bilinear fermion operators. 
It is not difficult to carry the calculation of one-particle matrix elements in p3| over to the present case of two 
components: 



N-l r 2ir r 2n j j 

^ / 2 c i(n-l)(u+ie)-i(q-l)(v-ie) /^-i^j (u)+4j (v) c ~»4v («)+»<M«) \ _ ^9) 



4tt 2 



Using the bosonic Wick theorem, the expectation value ( ) = exp [—Wo-] on the r.h.s. can be evaluated. At zero 
temperature, the function W a is given by 
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W„ = W„(u, v) = ^2 X! — ^ nw ~ amv cosm ( u + v )} _ costo(u - v)} . (30) 
This quantity is independent of component label <j, as expected. 

Comparing ( |30| ) with (39) in p3, it is seen that the effective coupling constants in the two-component case are 



1 2 1 2 

n ^ ] C^mi/j Tin 7) ^ ] Imu- 



W is a real and an even function of its arguments leading to the symmetries 



(31) 



(32) 



and to the condition n + q — 2m, m = 0, 1, 2, ... . 

For the interaction model IM1, one of the integrations in ( |30| ) can be performed giving the closed expression for the 
matrix elements of each component: 



M(m,p) = (c+_ p c m+p ) = -S pfi - — I ds 



sin((m + 1/2 - N)s) 



2sin(s/2) 

x exp[— 271(1 — cos(s))] 7p(2ai(l — cos(s))). 
Due to the factor {sin(...)}, the following symmetries hold: 



Similarily, IM2 leads to: 



M(m,p) = -6 pfi 



^2N—l—m—pcrC2N—l—m+p<T) — $p,Q ( c m—p o&m+p a) ■ 

^ dt cos(pt) f n ds \ sin((m + 1/2- N)s) 

2sin(s/2)) 



2tt [1 + Z a - cos(t)] a <J J ^ 2ir 



with decay parameters 



1 + Z 7 — cos(s) 
Zj = cosh(r 7 ) — 1, Z a = cosh(r Q /2) — 1. 



[(1 + Z a - cos(i - s))(l + Z a - cos(t + s))] 



So/2 



(33) 



(34) 



(35) 



(36) 



IV. NUMERICAL RESULTS 



The main results of the paper are the formulae ( |33| ) and ([35|) for the one-particle matrix elements. They are identical 
in form to those in |p4j| . depend, however, differently on the coupling constants. This leads to very different physical 
predictions, which are presented in a number of figures for fermion numbers 27V = 14 + 14. 
Using ( |l7|) and (|2p|), it is found that the main coupling parameters a m 



Vim) 



y/l + 2V{m) yjl -2V(m) 



(37) 



are non-positive and even functions of the interactions V(m): Irrespective of the sign of the interaction between the 
two components, the effective interaction in each component of the Fermi gas is attractive. 

In the case of IM1, only a\ is needed as input parameter in the calculation of the matrix elements. | V(l)| is obtained 
via (|3^) and all other quantities such as (i„ and 71 can be calculated from (|||) and (|3l|). 

We start with the discussion of the occupation probabilities P(m) = M(m,p = 0) of oscillator states as shown in Fig. 
1. It is seen that interactions smooth out the Fermi edge at mp = N — 1, but still leave a gap (not an energy gap!) 
at wiF- 
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FIG. 1. Occupation probabilities P of oscillator states m — 1 (m = 1, 2, ...) for an interacting two-component Fermi gas of 
2N = 14 + 14 atoms in a one-dimensional harmonic trap at zero temperature. Interaction model 1 with Qi = — 1 has been 
used. 

Fig. 2 displays the off-diagonal matrix elements for p — 1 . They are significant near the Fermi edge mp = N — 1 and 
cannot be neglected. Their values increase further with increasing coupling strength. 
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FIG. 2. Off-diagonal matrix elements M) versus oscillator state m — 1 (m = 1,2,...) for an interacting two-component Fermi 
gas of 2N = 14 + 14 atoms in a one-dimensional harmonic trap at zero temperature. Interaction model 1 with ai = —1 has 
been used. 



We also present results for the particle density and the momentum density. Both are expected to show Friedel 
oscillations [[ll) as noted in @,|4 1 . In agreement with , the effective intra-component interaction, which is always 
attractive, suppresses the Friedel oscillations in the particle density 



00 m 



m—0 p— — 7n 



(38) 



as is seen in Fig. 3. In (38), tp m (z) is the oscillator state |m) in position representation. 



G 



n(z)l 




FIG. 3. Dimensionless particle density n(z)£ (I is the oscillator length) versus dimensionless distance z from the center of 
the one-dimensional harmonic trap for 2N = 14 + 14 atoms in the two-component Fermi gas at zero temperature. Broken 
curve shows unperturbed Friedel oscillations. Smooth curve refers to the interacting case with ai = —1. Interaction model 1 
has been used. 

Conversely, the Friedel oscillations in the momentum density 



m—0 p— — 7n 



(39) 



are enhanced [Q. This is displayed in Fig. 4 for strong coupling [a\ = —10). We have chosen the oscillator length 
I = \fhj ' (jriA LOi) as unit of length, rendering n(z) and z as well as p(k) and k dimensionless. 



P(k)/1 




FIG. 4. Dimensionless momentum density p(k)/l {I is the oscillator length) versus dimensionless momentum k for 
2N = 14+14 atoms of a two-component Fermi gas in a one-dimensional harmonic at zero temperature. Broken curve 
shows unperturbed Friedel oscillations. Thick curve refers to the interacting case with «i = —10. Interaction model 1 has been 
used. 
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In the case of IM2, some modifications occur. We again set d% = —I. We need a m and 7m (m = 0, 1,2...) for the 
evaluation of This requires the knowledge of the decay constants r a and r 7 (cf. For convenience, we 

set r a = r 7 = r and estimate r by the following argument: The minimum wave number increment in the trap is 
Afc w 1/Lp oc l/ViV, where Lp — \/2N — 1 is the half- width of the classically allowed region at the Fermi energy. 
We, therefore, set r sa 1/yN or roughly r = 0.3 for the present case N = 14. This gives da = —1.16 for d\ — —1 and 
7o = 1.19. 

Fig. 5 shows the occupation probabilities of oscillator states for IM2. It is seen that they are more smoothly distributed 
than in the case of IM1, but still leave a gap at the Fermi edge. 

P(m-l) 
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FIG. 5. Occupation probabilities P of oscillator states m — 1 (m = 1, 2, ...) for an interacting two-component Fermi gas of 
2N = 14 + 14 atoms in a one-dimensional harmonic trap at zero temperature. Interaction model 2 with Qo = —1.16 has been 
used. 



Finally, we show the momentum density for IM2 in Fig. 6. The Friedel oscillations are still recognizable for small 
momenta, but strongly suppressed for momenta approaching kp = \/2N — 1. 



P(k)/1 




FIG. 6. Dimensionless momentum density p(k)/l [l is the oscillator length) versus dimensionless momentum k for 
2N = 14+14 atoms of a two-component Fermi gas in a one-dimensional harmonic at zero temperature. Broken curve 
shows unperturbed Friedel oscillations. Thick curve refers to the interacting case with ag = —1.16. Interaction model 2 has 
been used. 
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The off-diagonal matrix elements are significantly smaller for IM2 than for IM1. Nevertheless, they cannot be ne- 
glected: By comparing (|3^) with (|39|), it is seen that particle and momentum density would coincide in such an 
approximation. 



V. FERMI EDGE 



Fig. 1 and also Fig. 5 do not show the gapless distribution of occupation probabilities near the Fermi edge mp = N — 1, 
which is characteristic of a Luttinger liquid, i.e., our system is not a Luttinger liquid. This cannot be expected because 
the system is finite. 

We can, however, get a glimpse at Luttinger liquid behaviour in a special limit, which also presupposes a large particle 
number N . 

First, we consider a very slow decay of the interaction modes V(m) in IM2, i.e., r a — > r 7 -C 1. The factor 



1 + Z 1 — cos(s) 



7o 



7o 



(40) 



in the large square brackets of the integrand in (]35|) then becomes sharply localized at s = 0. 
We now calculate the occupation probability P(Ak n ) 



P(Ak n ) 



(41) 



near the Fermi edge and for N 3> 1. P(Ak n ) becomes a quasi continuous function of the wave number deviation 
Ak n = k n — kp = n/Lp — ► Ak, provided \n\ <C min(iV, l/r 7 ) is fulfilled. Using (40) in ( |35| ) we obtain 



P(Afc) = - - 



/ 1 3 

3F2 7o,7T,l;l,-;- 



r 7 L F 



Ak 



(42) 



in terms of a generalized hypergeometric function. It is seen that P(Ak) depends linearly on the wave number 
deviation in a small region near the Fermi edge. This can be compared with the Luttinger liquid prediction (cf. e.g., 



P LL (Afc) = i-sgn(Afc)C|A£f. 
C is a constant and the exponent [3 depends on the Luttinger liquid coupling strength "fpp according to 

(3 = 2j LL for 7 LL < i, 



(43) 



(44) 



and 



(3=1 for 7 LL > 



1 



(45) 



We conclude that the above limit of our model agrees with the case 7^^ > 1/2 of the Luttinger liquid. 

VI. DISCUSSION AND SUMMARY 

For the interaction to become significant in the quantities calculated, its strength V(l) should be as large as |U(1)| « 
0.5. We demonstrate that this condition is within experimental reach. 

To this order, we consider the dipole-dipole interaction p"5[ . It is marginally long ranged and thus favors forward 
scattering. In J3jj], it is shown that the inter-component interaction between longitudinally aligned dipoles reduces 
exactly to the effective one-dimensional potential 



Vm{k) 



2tt 



2af 



Ei 



2al 
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in momentum space. Here, at is the inverse of the transverse oscillator length, [i the magnetic dipole moment, and 
Ei denotes the exponential integral. 

Using this equation in the exact formula (A. 13) in p^ij ], V(l) for N = 14 is found to be 



V(l) = 0.8 



2 3/2 1/2 \ -, 
ml UJ e ' \ 1 

2^fF 2 J F' 



The quantity F denotes the filling factor F = Nui(,/u)f For example in 53 Cr, V(l) becomes of the required magnitude 
provided F is very small, i.e., the trap is highly anisotropic. 

In summary, the bosonization method has been used to construct a theory for a two component gas of spin polarized 
fcrmions in a one-dimensional harmonic potential with forward scattering between the two components. Asymptotic 
results with respect to the fermion number N were obtained for the one-particle matrix elements and used to discuss 
occupation probabilities for oscillator states, off-diagonal matrix elements, and distribution functions for particles 
and momenta in the harmonic trap. All these quantities can be significantly affected by the attractive interaction 
generated within each component. Specifically, the Friedel oscillations in the particle density are suppressed, while 
they survive in the momentum density. 

It has to be seen, whether the predicted Friedel oscillations can be observed experimentally. The amplitudes of the 
Friedel oscillations scale as 1 /N pij , hence Friedel oscillations are unobservable in a macroscopic bounded Fermi sea. 
Small particle numbers pose, however, severe detection problems. A conceivable experimental method to observe 



Friedel oscillations for atom numbers of the order of 100 is indicated in |34|. The method proposes microfabrication 
techniques to produce arrays of microtraps. 

On the other hand, the asymptotic bosonization method requires particle numbers, which are not too small. This 
is due to the presence of the anomalous vacuum, which couples to the real particles. For instance, the sum rule 
J2 n P(n) = N gives a somewhat larger value than the number N of real particles when (^3|) or (|55|) are used. The 
excess AA^ > grows with coupling strength and decreasing particle number. For N = 14 and very strong coupling 
a x = -10, AN is about 8 • 10~ 3 , for a x = -1, AN is less than 10 -10 . 

The atom number 2N = 14+14 and the coupling values employed here are appropriate to give visible Friedel 
oscillations and reliable results of the bosonization method. 

Acknowledgements: The authors thank S. N. Artemenko, F. Gleisberg, and W. P. Schlcich for valuable discussions 
and the Deutsche Forschungsgemeinschaft for financial support. 



VII. APPENDIX 

In this Appendix, we extend the bosonization procedure in p0| to the case of two components. 

Instead of d p ± operators, which are needed for the diagonalization of the interacting Hamiltonian, the following set 
of operators play the role of the b and b + operators (n > 1) in |40| |: 

l>na = [<%+ + Ka = -J= K+ + O^n-]- (A.l) 



They are canonical conjugates. Evidently: 



bL = -^PAn)- (A.2) 

Jn 



Then the two relations hold 



feU cX a ,c lal \ = Sa,^ -^= {ct +na ci a - c^q_„ ct ), (A.3) 

[bncr, cX a iho'\ = 5<j,a' — T= (Cfc_ no .C i(T - c£ ff Cj +no .). 



Following the arguments in KB], the associated Bose fields for a — a' are 



r (v) =- 1 V T e mv b na E-iV — e mv (d n+ + ad n .) f 4>+{v). (A.4) 
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